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Abstract 

We calculate the finite volume mass gap M(L) at 3-loop level in the non-linear 
O(n) a- model in two dimensions in small volumes. By applying the Monte 
Carlo measurements of the running coupling g 2 (L) = 2nM(L)L/(n - 1) by 
Liischer, Weisz and Wolff |I| measured in units of the physical mass gap m, 
the result is used to determine m in units of the A-parameter in the 0(3) 
and 0(4) models. Our determinations show good agreement with those by 
Hasenfratz, Maggiore and Niedermayer [j| in both models. We note that 
this manuscript has been revised in our paper hep-lat /98f 0025 by using the 
corrected four- loop (5- function on the lattice. 



1 Introduction 



The non-linear O(n) cx-model which describes an n- component spin field has thus 
far found many applications in theoretical physics. In condensed matter physics it 
has been investigated to study, e.g. ferromagnets. In elementary particle physics, 
on the other hand, the model in two dimensions shares with QCD many common 
properties; on the quantum level it is like QCD renormalizable |3J, asymptotically 
free according to perturbation theory P, || |7| and thought to have a mass gap 
m. In addition, the 0(3) model has instanton solutions. Further, the a-model 
has an additional advantage: because of low dimensionality it is simple to simulate 
numerically. Therefore, in many cases the model can be used for tests of new ideas 
in lattice theories. 

One of the most interesting properties of the a-model is the existence of an in- 



finite set of conserved non-local quantum charges |8], || [TO], plj . From these, one 
can exclude particle production and also derive the "factorization equations" di- 
rectly. Using these properties, one can then determine the S-matrix exactly up to 
CDD-ambiguity [12], provided the particle spectrum of the model is known Jl3, 14] . 
Furthermore, the CDD-ambiguity can also be restricted with help of the assumption 
that there are no bound states. The absence of the bound states was shown in the 
1 /^-expansion to the order of 1/n [13|. Further, the exact S"-matrix constructed in 



this way was confirmed by the same method to the second order in 1/n |~5 



Since there is only one scale in the model, one can, in principle, express all 
physical quantities with help of this scale. In particular, the mass gap of the theory 
is proportional to the A-parameter. This situation is similar to what happens in 
QCD with massless quarks where the A-parameter sets the scale for scale violations, 
for example, in deep inelastic scattering. 

Using the exact S'-matrix, a few years ago Hasenfratz, Maggiore and Niedermayer 
succeeded in an analytic determination of the ratio m/A in the 0(n) a- model |2], |3|. 
In this derivation, however, various assumptions including the thermodynamical 
Bethe ansatz are made, which are plausible, for which, however, the rigorous proofs 
are still lacking. An independent calculation of m/A is therefore desirable. 

In the literature 0, Liischer, Weisz and Wolff estimated the ratio m/A by mea- 
suring; the coupling g 2 (L) = 2nM(L)L/(n - 1) running; with the volume L of the 
system in units of the physical mass gap m[=M(oo)] by means of Monte Carlo 
simulations and applying the measurements to the perturbative expression of A- 
parameter in this scheme. Although their determination of m/A where A was eval- 
uated up to three loops in perturbation theory reproduces that of Hasenfratz et. al 
quite well, the agreement of two results was still not completely undisputable. In 
order to determine m/A more precisely, we decided to evaluate the A-parameter to 
one loop order higher. This requires the calculation of the finite volume mass gap 
M(L) at 3-loop level. 

Our work is also interesting in another respect. It is generally accepted that the 
0(n) a-model is asymptotically free. This assumption was, however, questioned by 
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Patrascioiu and Seiler [|T6| , [I~7| , Since in the determination of mj A by Hasenfratz 
et al. the validity of the asymptotic freedom was required, the confirmation of their 
result through our independent determination will, at the same time, support the 
correctness of asymptotic freedom. 

The work is arranged as follows. In the next chapter, we present our 3-loop 
calculations of the finite volume mass gap M(V). By computing the spin-spin cor- 
relation function, we first evaluate the mass gap on a lattice. We then take the 
continuum limit and finally convert it to the MS-scheme. In chapter ^, we apply 
the result to the determination of mass gap m in units of the A-parameter in the 
0(3) and 0(4) models. Two approaches will be considered. First, we introduce 
the mass shift 5q = \M(L) — m}/m which is known for L — > oo and match this 
with the perturbative M(L) |H| 2D , TL \. As a second approach, we determine m 



by applying the Monte Carlo measurements of the running coupling g 2 {L) to the 
perturbative A-parameter JTJ. Finally, the Feynman diagrams and some technical 
details are presented in the appendices. 



2 Calculation of mass gap at 3-loop level 

2.1 Mass gap on a lattice 

We consider an n-component spin field q l (x) (i — 1, • • • , n) with unit length q(x) 2 = 
1 on a two-dimensional finite lattice^ 

A = lx E Z 2 | - T < x < T; x 1 = x x + lJ, (2.1) 

where T and L are fixed integers with T,L > 2. We work with periodic boundary 
conditions in space direction and, for reasons which will become clear below, free 
boundary conditions in time direction. 
The action is given by 

n L r T ^ T 1 

5=rrE E d q(x) ■ d q(x) + £ d iq (x) ■ d x q{x) , (2.2) 

A J0 Xl= i l. Xo =-T x =-T * 

where /o denotes the bare coupling constant. The forward and backward lattice 
derivatives are defined by 

d^qix) = q(x + jx) - q(x) 
d *M x "> = - q(x - fl). 

fi is the unit vector to the positive /^-direction (/x = 0, 1). The expectation values of 
the observables can be calculated by the formula 



(°) = \\ II [d n q(x)5(q(x) 2 - l)]Oe" 5 , (2.3) 



1 We set the lattice constant a to 1. 
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where the partition function Z is such that (1) = 1. 

Concerning the energy spectrum, there is a following possibility for the calcula- 
tion of the mass gap M(L) in finite volume for fo — > 0. We calculate the expectation 
value of the spatially averaged spin-spin correlation function in the limit T — > oo 
(r>0): 

C(r) = hm -L £ £ . g^U^. (2.4) 

If the bare coupling fo is fixed, the correlation function C(t) converges to the vac- 
uum expectation value of the 2-point function for spin field operators at large T 
independently of whether one chooses free or periodic boundary conditions in the 
time direction. In perturbation theory, we are, however, first expanding in powers of 
the coupling f and then let T go to infinity. Therefore, we can not, in general, have 
such a converging behavior, and one requires a proper choice of boundary conditions 
to get the desired vacuum expectation value. 

In free boundary conditions, we are projecting on states which are O(n) invariant 
at large times. The energies e(/o) in these states, on the other hand, have the 
property 

Bm{e(/o)-e (/o)}>0, (2.5) 

where £o(/o) denotes the ground state energy. Except for the ground state, all these 
states are therefore exponentially suppressed at large T and we arrive at the desired 
vacuum expectation value even though we are first expanding in powers of fo [0. 

After having taken T to infinity so that C (r) been reduced to the vacuum expec- 
tation value of the product of two spin field operators, only the vector intermediate 
states contribute. The energy £i(/o) of the first excited state of this system has the 
property 

£i(/o) - e„(/o) = O(f ) for f - 0. (2.6) 

For the energies of the other higher excited states, eq.( |2.5| ) is valid and their contri- 
bution to the spin correlation function hence vanishes exponentially at large r. 

From this consideration, the mass gap M(L) which is defined by the l.h.s of 
eq.( |2.6| ) can be calculated by 

M(L) = -\im 1 -^-\nC(r). (2.7) 

If one expands C(r) in perturbation theory for fo — > 0, the mass gap has the general 
form 

1 oo 

M{L) = —Y^K^ ■ (2-8) 

can thus be determined by calculating the 2-point function |eq.( |2.4j )1 with the 
mentioned boundary conditions in perturbation theory. 
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The calculations up to the third order in fo were already done by Luscher and 



Weisz |]22| . Their results read 



A (D 

A (2) 
A (3) 
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fc 
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0.11141943681128(1) 
0.001927404 148 69(1) 



(2^FV ln ^- 7 



(2.9) 
(2.10) 
(2.11) 



(2.12) 
(2.13) 

(2.14) 
(2.15) 



(2.16) 
(2.17) 



(2.18) 
(2.19) 

(2.20) 



7 denotes the Euler constant (7 = 0.577216 • • ■). 

In order to determine the mass gap in the fourth order, we first evaluate A^ 4 ) on 
finite lattices and then extrapolate to continuum limit, i.e., L/a — > 00. The 2- and 
3-loop diagrams contributing to A^ 4 ^ are illustrated in figures (§) and (g) which are 
presented in appendix A due to the very large number of the diagrams. In earring out 
Wick's contractions from them, we used the symbolic language MATHEMATICA. 

We calculated the generated terms numerically. In the numerical computations, 
we take T and r large, but finite. The corrections to the limit T, r — > 00 are here 
exponentially suppressed by the order of e~~^ 2T and e~^^ T ~ T \ We therefore achieve 
the best approxiamation to the limit T, r — > 00 by keeping T ~ 3r. 
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The numerical work, however, turned out to be very complicated due to the 
problems caused by rounding errors and running time which increases with L very 
quickly (in the order of L 6 ). In addition, problems appear from the fact that on the 
3-loop level on which we are working there are extremely many terms to treat. One 
needs therefore a very efficient computer program which requires, among others, the 
free propagator running fast without significant loss of precision. In this way, we 
succeeded in evaluating all diagrams up to L = 20 with good enough precision and 
reasonable running time (for L = 20 we needed the CPU time of around 7 days on 
the IBM RISC 6000/32H), so that we could extrapolate to the continuum limit to 
get our desired results. 

With regard to the extrapolation to the continuum limit, we note that in this 
limit A^ 4 ) has the general form 

A< 4 > = A[, 4) + Ai 4) (lnL) + A^QnL) 2 + A^QnL) 3 (2.21) 

up to terms of the order L~ 2 (lnL) 3 where the coefficients Ao 4 -*, • • • , Ag 4 "* are inde- 
pendent of L. All of these coefficients except Aq 4 "* can be derived with help of the 
renormalization group equation and the results from the calculations of the mass 
gap at lower orders: 



(27r)n 4 
(n- l) 3 

(27r) 3 n 4 I'" 1 "\~~ y/2 



{/ii - 2(hi-^= - 7 ) + 3(2tt) 2 (A; 1 - A; 2 )| 
- 2/ii + h 2 + 4^1n -y= - 7) + vr + 3(2tt) 2 (A;o - fci)} 



(2vr) 3 n 4 I 1 ' z ' V y/2 
n — 1 f , , /, 7r 



' - A. - 2(ln^ - 7) — 3(2,) 2 *o} (2.22) 



(27r)% 4 l V V2 

A?» = -^?(ln4-7 



(27r) 3 n 4 I 
(n- l) 3 



{ 9 ( ln * 7 ) + I(3, + 5)} 



(2vr) 3 n 4 I V v 7 ^ 
(n-1) 



{3(31n^ + .-3 7 + l) + 2} 



(27T)% 4 I V y/2 

+ ^{ 1 + l{ 2 ^7=2 + n - 2l + 1 )} (2 ' 23) 

(4) (n - l) 4 _ 3(n-l) 3 3(n-l) 2 _ n-1 

3 (27r) 3 n 4 (2vr) 3 n 4 (27r) 3 n 4 (2vr) 3 n 4 1 ' ' 

In order to determine the unknown coefficient Aq 4 -*, we decompose it into n- 
independent constants: 

A o 4) = ^I^{ s ° + ( n ~ X ) s i + ( n " l ? s * + (n- l) 3 s 3 }. (2.25) 
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The coefficients sq, • • • , S3 can then be determined by inserting A^ 4 ) which was eval- 
uated for different L's in eq. ( [2.21 ) and extrapolating to continuum limit. 

The extrapolation of Sfc (A; = 0, • - • , 3) to L — 00 can be done very efficiently 
with help of the procedure by Liischer and Weisz described in detail in ||23|| . We 
note that in our case the constants Sk have corrections of the form 



s°(L) = Sfc + ^{a p + 6 p lnL + c p (lnL) 2 + rf p (lnL) 3 }/L 2p . (2.26) 



Through the extrapolation in the region 5 < L < 20, we obtain 



s 
s 3 



0.03954(1) 
0.02903(1) 
0.000756(1) 
-0.000649(1) 



(2.27) 
(2.28) 
(2.29) 
(2.30) 



At this point, we would like to emphasize that for all four constants the limit 
L — > 00 exists. This is, of course, only the case if A^ 4 ) calculated by us contains the 
terms diverging logarithmically with L, so that they cancel exactly with the other 
divergent factors of eq.( |2.21|) in a non-trivial way. This correct logarithmic behavior 
of AW is a strong consistency check on our evaluations of the diagrams. 



2.2 Conversion of the mass gap to the MS-scheme 

The /^-function in the MS-scheme of dimensional regularization is defined by 



(2.31) 



where / denotes the renormalized coupling and \x the normalization mass. The first 
two coefficients in (|2.31| ) are scheme independent: 



bi 
b 2 



71-2 

n-2 
(27m) 2 



(2.32) 
(2.33) 



The remaining coefficients are known to four loops |P4] , p5| , |26| , p7 



b, 



n-2 



(27rn) 3 

n-2 

(2t™) 4 



(n-2)- + l 



n-2) 2 — + (n-2) Ul + u 2 



(2.34) 
(2.35) 
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with 



Ul 
"2 



§(«+!) 

-i(3,-i; 



(2.36) 
(2.37) 



where the constant v is given by v ~ 1.2020569. 
The corresponding A-parameter is 

Ams = Kbifr b2/bi e- 1/M) ■ X(f), 

rf 



Kf) 



cxp 



dx 



1 1 

+ 



o \(3(x) bix 2 b\ 



x 



(2.38) 
(2.39) 



The function A(/) is well-defined at / = and can therefore be expanded as a power 
series in /. Up to 4-loop, we find 



A.w., M&i/)- 62/6? e- 1/(6l/) U + ^|^/ 



b\ 

b\ - b\b\ + 2blb 2 b 3 - 2b x b% + gjj - frfo 4 
26? 



f + o(f 



On the other hand, the (3 function on the lattice is defined byQ 



fcfo) 



da 



-fa E bJo ■ 



f v=\ 

The coefficients here are also known to four loops f28| , |29"| , |30| : 



6i 
63 



bi 
b 2 

n-2 
(27m) 3 
n-2 
(2vrn) 4 



[(n - 2)/ii + h 2 ] 
(n - 2) 2 t! + (n - 2)£ 2 + < 3 



with 



/i! = -0.088766484(1) 

/i 2 = l + 7r/2-57r 2 /24 

*i = -1.015(1) 

t 2 = -5.44(1) 

t 3 = -9.093756(4) 



(2.40) 



(2.41) 



(2.42) 
(2.43) 

(2.44) 
(2.45) 



(2.46) 
(2.47) 
(2.48) 
(2.49) 
(2.50) 



2 In this section, we introduce the lattice constant a. 
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where the difficult 4-loop computation was performed recently by Caracciolo and 
Pelissetto [[3(J . 

The A-parameter on the lattice can be obtained from Ams of eg. Q2.38] ) 
through the replacements p — > a -1 , / — > fo and f3(x) — > The two A-parameters 

are related by the formula: 



A, 



A 



MS 



7T 



exp - < In 7 



7T 



n-2 



(2.51) 



From this, it follows the relation between the coupling constants in both schemes. 
Up to order C(/ 4 ), we find 
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X (3) 



, oo 



n-2 
Airn 
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4n 
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(2.52) 

(2.53) 
(2.54) 

(2.55) 



Finally, we use this relation to express the mass gap [eqs. (|2.9|) - (|2.11|) and 
( |2.21| )] in terms of the renormalized coupling in MS-scheme. After a lengthy, but 
straightforward calculation we get 



M(L) 
E (i) 

-(2) 



n 
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n-2 
4irn 



/ i + E ^ v) f u 

, P 2 L 2 
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2(n 
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'3 (47m) 



(47m 



(47m) 2 
n-2 



-K 



(1) 



(47m) 



(47m) 



(2.56) 

(2.57) 
(2.58) 



(2.59) 



The expressions for the coefficients pi, - ■ ■ , Pq of k^ 3 ' are rather long; we therefore 
write them in appendix B. Here, the coefficients and are the result of Liischer 
and Weisz |22|]. By using dimensional regularization, was first calculated by 
Floratos and Petcher pT] with less numerical precision. Our result for tffl is the 
extension of their 1- and 2-loop computations to 3-loop. 
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With dimensional regularization, k^ u ' = for n = 2 since we formally have a free 
theory in this ^ and k^' already satisfy this condition. We checked that 

k' 3 ' is also zero within the errors, which shows a further check on the correctness 
of our final results in the lattice calculations, s , si, s 2 and s 3 [eqs.( p.27|) -( p.30| )1. 
This is, however, not only a check on our computations, but also that on the 4-loop 
coefficient of the (3 function on the lattice 64 which was introduced in our calculations 
by the conversion of the result on the lattice to the MS-scheme. 



2.3 Expansion in z 

In view of the evaluation of the finite volume mass gap M(L) in the infinite volume 
limit, we introduce the variable 



z = M(L)L 



(2.60) 



and express the ratio M(L) /Ajj^ in terms of this dimensionless and renormalization 
group invariant parameter in perturbation theory for z — > 0, as suggested by Liischer 
19fl . Here, Aj^g is related with Ams of eq. fl2.38Q through 



A 



MS 



A MS exp 



(ln47r — 7) 



(2.61) 



For that purpose, we invert at first the /-expansion of eq.( 2.56 ) in the z- 
expansion: 

f = diz + d 2 z 2 + d 3 z 3 + d 4 z 4 + 0(z 5 ), 



where 
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The insertion of this series in eq. fl2.38 ), together with eq. fl2.6lQ , yields finally 
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(2.68) 
(2.69) 
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The expressions for the constants a%, ■ • ■ , 0:4 are given in appendix B. We note that 
the coefficients a u in eq.( [2.67|) are independent of /i and L since M(L), and z 
are renormalization group invariants. In the above result, the leading term and the 
first coefficient a\ were calculated in ref. 
3-loop correction to their calculations. 



221, and our determination of a 2 is the 



2.4 Definition of running coupling g 2 {L) 



Observing that the mass gap M(L) in finite volume |eq.( [2.56D l is proportional to / 
in the leading order, Luscher, Weisz and Wolff IJ defined a coupling running with 
/- by 

g 2 (L) = 2nM(L)L/(n - 1). 



We can then expand it, for p = 1/L, in the coupling of the MS-scheme: 

f(L) = f + c 1 f 2 + c 2 f + c 3 f 4 + 0(f), 



where 
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For the coefficients of /^-function 
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we find 
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(2t™) 3 
1 n-2 



4 (27m) 



n - 2) 3 X i + (n - 2) 2 X2 + (n - 2) X s + X4 



(2.70) 
(2.71) 

(2.72) 
(2.73) 



(2.74) 



(2.75) 



(2.76) 
(2.77) 

(2.78) 
(2.79) 



where we have used the /3-function in MS-scheme |eqs. (|2.3l]) -( 2.35 )l. The expressions 
for the constants \ii ' ' ' > X4 are listed in appendix B. 
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The corresponding A-parameter is given by 



A(# ) = exp 



hhfy^e' 1 '^ -\{g% (2.80) 



/ dx( 

Jo \8(: 



(3(x) hx 2 bjx 



(2.81) 



From the relation between two couplings in eq. (|2.71| ), we can derive the ratio of 
A-parameters in both schemes: 



A FV 



Ams 2^ 

3 Results and Discussions 



(2.82) 



In this chapter, we determine the mass gap m in units of the A-parameter by ap- 
plying our results in the last chapter. For this, the following two approaches will be 
considered: 

• Determination of the ratio m/Ajg-g by matching the behavior of the finite 
volume mass gap M(L) evaluated at small L with that of the mass shift 
5q = [M(L) — m]/m known at large L 

• Determination of the ratio Apy/m by means of running coupling g 2 {L) 

We discuss these two methods in detail and give the results. 

After that, it is interesting to compare them with the analytic determination 
performed by Hasenfratz, Maggiore and Niedermayer a few years ago [0, [3]. Their 
result in the general O(n) model has the following simple form: 

/8\V(»-2) 1 

m= \- e ) r[i + i/( B -2)] A ** (3A) 

In our work, we will need the values in the 0(3) and 0(4) models: 

m = for 0(3) (3-2) 

/32\ V 2 

m = (-J A^ for 0(4) (3.3) 



By using the equations ( |2.82j ) and ( J2.61 ), these formulas can be rewritten in 

with the running coupling g 2 {L): 

for 0(3) (3.4) 



terms of Apy defined in eq. ( |2.80|) in connection with the running coupling g 2 (L): 

Apy e^ +1 



m 327r 

— = —\[^e~< for 0(4) (3.5) 
m 16 V 2tt w ( j 
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3.1 Determination of m/Aj^g by applying 5q 

Our aim is to determine the mass gap m in infinite volume by applying our result 
of the mass gap M(L) in finite volume [eg. ([2.671) 1 . F° r this purpose, we introduce 
two quantities C (z) and c by 

M(L) = C (z)A m (3-6) 
m = c o% ( 3 - 7 ) 

where m is obtained from M(L) by m = lim^oo M(L). 

The problem is then the determination of Cq which may be done by attempting 
to extrapolate C (z) to the infinite volume limit, i.e. to z — > oo. Inspired by the 
successful application of the method in the O(n) model at n = oo where one can 
solve the given problem exactly, Luscher |19| tried, with his 1-loop calculation of 
Co(z), to estimate Co in the 0(3) model as the value of Cq(z) at some intermediate 
region of z (at around 3 < z < 4) with the hope that in this region the perturbative 
Co (z) would be a good approximation for the infinite volume. Here, we extend his 
1-loop calculation to up to 3-loop which is plotted in fig. |] as a function of z in 
the 0(3) model. The lowest curve is the 1-loop result by Luscher while the middle 



one comes from the 2- loop computation by Floratos and Petcher [pi]] . Finally, the 
uppermost curve shows our 3-loop result. The expected behavior of the curves, 
supported by the study in the 0(oo) model, is that they rapidly decrease at small 
z and then quickly become flat. From the figure we, however, see that it is difficult 
to make a clear estimation for Co(oo) in this way. 



For a better estimation of cq, we make use of the mass shift defined by 

M (L) - m 



m 



(3.8) 



as done in ref. [plfl . 5q was first introduced by Luscher pOj and can be calculated 



at large £ = mL exactly up to an exponentially small correction term if the elastic 
forward scattering amplitude is known.0 In the 0(n) a- model, the exact scattering 



matrix was determined by brothers Zamolodchikov |13|, [TJJ]. In particular, the for- 
ward scattering amplitude is known and hence one can calculate 5q in this model at 
large C- 

By using the definition of 5q together with eq. (p.7|), we obtain the following 
relation between Cq(z) and c$: 

C = (l + 5 )c . (3.9) 
From this, a formula expressing Co as a function of ( follows: 

co(C) = [l + 5o(C)]" 1 CoKC)]- (3.10) 



3 £ is related with z by z = (1 + <5q)C so that z ~ £ at large £. 



12 




Figure 1: Co as a function of z in the 0(3) model. The lowest curve contains 
only the 1-loop contribution, the middle one the contributions up to 2-loop and the 
uppermost curve up to 3-loop. 

Actually, Co is independent of (. Therefore, it must be valid for all (: Cq[z(C)) ~ 
1 + Sq(Q if Co and So are known exactly. However, we know the two quantities only 
approximately; Co for small z and So, on the other hand, for large (. Although in this 
approximate relation Co is in general dependent on (, there may exist the possibility 
that one finds somewhere an overlapping region where both approximations are good 
and thus Co becomes independent of (. 

This idea works very well in the model for n = oo, as shown in ref. jnj]. We 
therefore apply it to the determination of Co in the 0(n) model, explicitly in the 
0(3) and 0(4) models. In the 0(3) model, Sq was calculated by Liischer pO| : 



So = 4vr r dte-^^±- 2 + 0(e~% (3.11) 

where k is a constant which is not smaller than y| and may become as large as 3. 
If we now insert this quantity together with Cq(z) at n = 3 in eq.( ^.10| ), we obtain 



a functional relation of Co in terms of ( which is illustrated in fig. 0. 

The 1-loop result displays a region where Co — 1.6. This value was the earlier 



estimation by Liischer by this method fTS, BO]. There, Liischer gave an optimistic 



error of 20 % which he apparently underestimated. The 2-loop result corrects this 
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Figure 2: cq as a function of ( in the 0(3) model. The lowest curve contains for 
Cq{z) only the 1-loop contribution, the middle one the contributions up to 2- loop 
and the uppermost curve up to 3-loop. 

by about 30 % to Cq ~ 2.15. This is in agreement with the estimation by Floratos 
and Petcher who give cq ~ 2.1 f2~H| . Finally, we see from our 3-loop result that 
the curve rises rather monotonically instead of showing a wider flat area. If this 
method should work well, we would expect that the flat area at the 3-loop becomes 
wider than at the 2-loop. We attribute the bad applicability of this method to the 
fact that the convergence radius of Cq(z) is apparently too small. Nevertheless, we 
conservatively estimate the 3-loop approximation to Cq by cq ~ 2.3. We note that it 
is very difficult to give a systematic error on our estimation. 

Our determination is to be compared with that by Hasenfratz et al. which is 
c = 2.94304- •• [see eq.( [3.2|) 1. The value of c — 2.3 at the 3-loop still deviates 
noticeably from that by Hasenfratz et al.. One does not, however, need to be 
disturbed because, from the convergence of the 1-, 2- and 3-loop approximations, 
one sees a definite indication that the results in higher orders would arrive at that 
of Hasenfratz et al.. 

Now, we investigate the 0(4) model and see how the picture changes in this 
model. The mass shift here has the form 

1 r°° ( 
5 = — dicoshte- Ccosh * 4- 

2lT J-oo I 



4 2tt ' 4 ~*~ 2ir > 



2nt a 



+2 + sl Bfi _ iL i + «/) 



e>( e - AC ), (3.i2) 
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Figure 3: Cq as a function of £ in the 0(4) model. The lowest curve contains for 
Cq(z) only the 1-loop contribution, the middle one the contributions up to 2- loop 
and the uppermost curve up to 3-loop. 



where B denotes the beta function and A is a constant whose value lies in the same 
region as k in eq. (|3.11| ). We again insert this quantity together with Cq(z) at n = 4 



in eq.( p.lO| ) to get the functional relation of Cq in terms of (. This relation is plotted 
in fig. H 

The 1-loop result shows a stable region in the neighborhood £ = 3, from which 
one may estimate Co — 1.5. If one, however, looks at the curves of the higher 
loops, they become steeper and rise rather monotonically, which makes it difficult 
to estimate the desired value. We again see, as in the 0(3) model, this method 
does not work very well. However, a remarkable information is contained in the 
estimations of Co from the curves; by going to the higher loops, they are tending to 
converge to the value by Hasenfratz et al. of c = 1.93577- • • [see eq. ([3.3|) 1- 

From our determinations of the mass gap m, we conclude that in both 0(3) and 
0(4) models our results support the correctness of those by Hasenfratz et al.. 

3.2 Determination of A_pv/ m by applying g 2 (L) 

In chapter 0, we introduced a coupling g 2 {L) running with L. There, the coefficients 
of the corresponding /3-function /3(g 2 ) were calculated up to 4- loop order |eqs. ([2.76|) - 
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771 L 




0019(3) 


0.5372 


0.0038(3) 


0.5747 


0.0063(4) 


0.6060 


0.0127(4) 


0.6553 


0.0211(5) 


0.6970 


0.0327(4) 


0.7383 


0.0422(4) 


0.7646 


0.0659(6) 


0.8166 


0.1317(10) 


0.9176 


0.2744(15) 


1.0595 


0.5557(13) 


1.2680 



Table 1: g 2 {L) in the 0(3) model as a function of L in units of m 



(|2.79|) 1. We make use of these perturbative coefficients to evaluate the A-parameter 
Apv in this scheme in perturbation theory. 

In order to obtain the perturbative expressions for Apy, we expand A(g 2 ) of 
eq. ([2.81|) in g 2 . Up to order 0(g 4 ), we find 



A (2) 

A 



FV 
(3) 



I(6 1 ^)-6 2 /6f e -l/(6 1 ^) 

L 



FV 



1 

I 1 



hg 



i + 



u 2 



hh-2 

g 



A 



(4) 
FV 



^b 1 g 2 y b ^e- 1 ^l [ l + 



b\ 9 



b\ 



b\b\ + 2blb 2 b 3 - 2b^% + b\b 2 - b% 



2b\ 



(3.13) 
(3.14) 



(3.15) 



where Apy, Ap V and Ap V represent the 2-, 3- and 4-loop approximation to Apy 



respectively. In the limit L 

converge to the exact 



where the coupling g (L) also goes to zero, they 



A 



FV 



2,-b 2 /bl e -l/(b 1 g 2 ) 



(3.16) 



with the rates g 2 (L), g 4 (L) and g 6 (L) respectively. 

For determination of these perturbative A-parameters, the values of the running 
coupling g 2 {L) are needed. In ref. 0, they were measured numerically by Liischer, 
Weisz and Wolff through a finite size technique for given L in units of the mass 
gap m in the 0(3) model. Furthermore, Hasenbusch made available to me his new, 
more precise data |31j which contain the measurements in the smaller region of 
the coupling. Their results are listed in table |l|. In deriving the numerical values 
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of the table, the extrapolation to the continuum limit was done by fitting with a 
polynomial in (a/L) 2 which was supposed by Symanzik [32 . 

If we insert the numbers of table [I] in eqs.( |3.13|) - (|3.15| ), we get A FV /m in 2-, 3- 
and 4-loop approximations. The results are listed in table || and also illustrated 
in fig. |j. In the figure, we also draw the analytic value of Hasenfratz et al. for 
comparison. We see in the figure the difference between the 3- and 4-loop results is 
very small. Nevertheless, our 4-loop result lies between the 2- and 3-loop ones, so 
that by consideration of the higher order for Apv the ratio Apyjm lies nearer to the 
analytic value. The remarkable fact, however, is that the 3- and 4-loop curves which, 
with the old data by Luscher et al. alone, seemed to fall down continuously again 
rise together to the compared value if one adds the new data for smaller coupling 
by Hasenbusch beginning from g 2 {L) = 0.6970. We can thus expect very well that 
in the limit g 2 (L) — > the perturbative A FV /m would converge to the result of 
Hasenfratz et al.. 

To give a further impression on the systematic errors, we consider an alternative 
definition of ra-loop approximations to Apv] instead of expanding X(g 2 ) in g 2 , we 
insert the perturbative coefficients of (3(x) in X(g 2 ) and integrate this exactly. If 
we apply the coefficients calculated up to 4-loop order, we obtain the following 
approximations to Apv at 2-, 3- and 4-loop levels which are a little modified from 
eqs.(FT3D-(gl5|): 



A (2) 



FV 



1 

t 



'i9 2 ! 



• exp 



-*>2/fe?„-l/(6i9 2 ) 



dx 



o 



+ 



1 



x 2 (b\ + b 2 x) b\x 2 b\x 



U2, 



A 



(3) 
FV 



-{b l g 2 )- b ^e- 1 /^ 
L 



■ exp 



dx 



x 2 (bi + b 2 x + b 3 x 2 ) b\_x 2 b\x 



A 



(4) 
FV 



L 



■ exp 



dx 



+ 



(3.17) 



(3.18) 



(3.19) 



.x 2 (6i + b 2 x + b 3 x 2 + b 4 x 3 ) ' hx 2 b\x, 

Insertion of the data of table [l] in these equations gives the numbers in table ^| 
which are also plotted in fig. [5]. Apart from the fact that the difference between 
the 3- and 4-loop results is a little larger, we have here almost the same picture as 
in fig. f|. Although the 4-loop result lies below the 3-loop one so that it deviates a 
little more from the value of Hasenftratz et al., there is no cause for alarm since both 
curves approach the analytic value together. From our study through the alternative 
definition of the perturbative Apv, we therefore conclude that our determination of 
Apv/ m , as in the first case, agrees with that by Hasenfratz et al. very well. 

In ref. [Q, Hasenbusch and Horgan measured the running coupling g 2 {L) as a 
function of L in units of the mass gap m in the 0(4) model by using the finite-size 
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q 2 (L) 


A ( ?l/m 

FY 1 


A { il/m 

Fvi 




0.5372 


0.0511(8) 


0.0467(6) 


0.0467(17) 


0.5747 


0.0512(7) 


0.0465(6) 


0.0466(17) 


0.6060 


0.0514(7) 


0.0464(6) 


0.0465(15) 


0.6553 


0.0518(6) 


0.0464(5) 


0.0464(15) 


0.6970 


0.0520(6) 


0.0462(5) 


0.0463(15) 


0.7383 


0.0525(6) 


0.0463(5) 


0.0463(13) 


0.7646 


0.0526(6) 


0.0462(5) 


0.0463(13) 


0.8166 


0.0532(6) 


0.0463(5) 


0.0463(11) 


0.9176 


0.0552(5) 


0.0471(5) 


0.0472(10) 


1.0595 


0.0574(4) 


0.0478(3) 


0.0479(8) 


1.2680 


0.0628(1) 


0.0502(2) 


0.0503(8) 



Table 2: 2-, 3- and 4-loop approximations to Apy/m in the 0(3) model 



0.08 
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1 0.8 
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0.6 



Figure 4: 2-, 3- and 4-loop approximations to Apy/m in the 0(3) model from 
table |2|. The triangular points show the 2-loop approximation to the A-parameter, 
the hexagonal points the 3-loop and the circular points the 4-loop approximation. 
The dotted line is the value of Hasenfratz et al. (Apy/m = 0.04815 • • ■). 
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Q 2 (L) 


r VI 


± 1^ y j 1 1 u 


A^l/m 

r V 1 


n 5^79 

\J .00 1 £ 


n n555f8'i 


U.U4: 1 4:1 U 1 


U.UtiUOI -L 1 I 


0.5747 


0.0559(7) 


0.0473(6) 


0.0466(17) 


0.6060 


0.0563(7) 


0.0473(6) 


0.0466(15) 


0.6553 


0.0572(6) 


0.0474(5) 


0.0465(15) 


0.6970 


0.0578(6) 


0.0474(5) 


0.0464(15) 


0.7383 


0.0586(6) 


0.0476(5) 


0.0465(23) 


0.7646 


0.0590(6) 


0.0476(5) 


0.0465(13) 


0.8166 


0.0601(6) 


0.0479(5) 


0.0466(11) 


0.9176 


0.0633(5) 


0.0492(5) 


0.0476(10) 


1.0595 


0.0671(4) 


0.0505(3) 


0.0484(8) 


1.2680 


0.0755(1) 


0.0544(2) 


0.0513(8) 



Table 3: Alternative 2-, 3- and 4-loop approximations to Kpyjm in the 0(3) model 



0.00 



0.07 



0.06 



\< 



0.05 



0.04 




Figure 5: 2-, 3- and 4-loop approximations to Kpyjm in the 0(3) model from 
table 0. The arrangement of the points and dotted line is same as in fig. f|. 
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mL 


S 2 W 


1/8 


0.863(2) 


1/4 


1.011(2) 


1/2 


1.228(2) 


1 


1.584(4) 


2 


2.309(10) 


4 


4.132(10) 



Table 4: g 2 {L) in the 0(4) model as a function of L in units of m 



scaling analysis of Liischer et al. pj. We make use of the measurements to determine 
the ratio Kpyjm in this model. Their data for the running coupling are given in 
table [§ If we apply the numbers in the table to the eqs.( [3.13| )-( [3.15| ), we obtain 
Kpyjm i n 2-, 3- and 4- loop approximations. They are written in table |5] and also 
illustrated in fig. |^. The figure shows very good agreement between our result 
and the one by Hasenfratz et al.. Already at g 2 {L) = 1.584, Kpyjm at 4-loop level 
is almost the same as the analytic value. 

As in the 0(3) model, by inserting the data of table f| in the eqs.( |3.17f )-( |3.19| ), 
we also determine Kpyjm by means of the modified definition of the perturbative 
Kpy. The results are written in table [5] and also plotted in fig. [7| Apart from the 
fact that the curves for the 3- and 4-loop approximations in the figure changed their 
places, we see here the same picture as in fig. || The agreement of two values at a 
little smaller g 2 (L) is irrelevant from the fact that we are at the end interested in 
the determination of Kpyjm at the limit g 2 {L) — > 0. 

As a whole, we conclude that not only in the 0(3) but also in the 0(4) model 
our determinations of Kpyjm agree with those by Hasenfratz et al. very well. 

Acknowledgement 

I thank Peter Weisz for suggesting this work and many useful discussions. He also 
read through this manuscript, which is greatly acknowledged as well. I express my 
thanks also to Martin Hasenbusch for having made available to me his unpublished 
numerical data. 



20 



?{L) 


kf v /m 


Kf v /m 




0.863 


0.0795(7) 


0.0722(7) 


0.0725(16) 


1.011 


0.0817(7) 


0.0729(6) 


0.0733(15) 


1.228 


0.0844(5) 


0.0734(6) 


0.0740(12) 


1.584 


0.0880(5) 


0.0733(6) 


0.0742(10) 


2.309 


0.0928(4) 


0.0701(4) 


0.0723(9) 


4.132 


0.0853(4) 


0.0479(4) 


0.0544(9) 



Table 5: 2-, 3- and 4-loop approximations to Apv/ m in the 0(4) model 



0.1 



0.08 - 



< 



0.06 



0.04 




Figure 6: 2-, 3- and 4-loop approximations to Apy/m in the 0(4) model from 
table [5[ The arrangement of the points is same as in fig. f|. The dotted line is the 
value of Hasenfratz et al. (Apv / m = 0.07321 ■ ■ •). 
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Ag^/m 


Afv/m 


0.863 


0.0830(7) 


0.0733(7) 


0.0726(16) 


1.011 


0.0859(7) 


0.0744(6) 


0.0734(15) 


1.228 


0.0898(5) 


0.0757(6) 


0.0743(12) 


1.584 


0.0952(5) 


0.0771(6) 


0.0749(10) 


2.309 


0.1035(4) 


0.0780(4) 


0.0741(9) 


4.132 


0.1024(4) 


0.0678(4) 


0.0609(9) 



Table 6: Alternative 2-, 3- and 4-loop approximations to Kpyjm in the 0(4) model 



0.1 - 



0.08 



0.06 




g 2 (L) 



Figure 7: 2-, 3- and 4-loop approximations to A^y/m in the 0(4) model from 
table p. The arrangement of the points and dotted line is same as in fig. |[ 
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A Feynman diagrams 



We list here the 2- and 3-loop Feynman diagrams contributing to the mass gap in 
fourth order. In the diagrams, the internal dot denotes the vertices while the cross 
means those coming from 'V 2 insertion." 



o 



(1) 



CG 

(2) 
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-X X- 
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o 



(10) 



e- 



(13) 
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a. 



(14) 



(12) 



(15) 



-X X- 



(16) 



Figure 8: 2- loop diagrams 
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Figure 9: 3-loop diagrams (continued) 

B Expressions of Coefficients 

In this appendix, we show the expressions of the coefficients defined in the main 
text. 



B.l Expressions of the coefficients appearing in eq. (|2.59| ) 
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B.2 Expressions of the coefficients appearing in eq.( |2.t)9| ) 
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B.3 Expressions of the coefficients appearing in eq.(|2779|) 
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